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An efficient method of accurately estimating delay spread in multi-path transmission applications is based on quantifying the effect 
of the actual delay spread on the shape of the correlation energy profile between the received signal (10) and the synchronization pattern of 
the receiver. Depending on the degree of the estimated delay spread (20), an appropriate demodulation technique (30 or 40) is selected for 
optimizing the receiver performance over a range of channel multi-path conditions. The invention is applicable to digital wireless mobile 
communication systems, as well as to any device that performs characterization or testing of a transmission system. 
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DELAY SPREAD ESTIMATION FOR MULTIPATH 
FADING CHANNELS 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims the benefit of Provisional Patent Application 
No. 60/1 1 1,703, filed December 10, 1998, entitled "Delay spread estimation technique 
in multi-path fading channels for mobile communication system," 

BACKGROUND OF THE INVENTION 

The present invention relates to digital wireless mobile communication 
5 systems employing a synchronization sequence, and more specifically, to methods and 
apparatus for estimating the delay spread in a time- varying multi-path fading channel. 

Modern wireless mobile receivers are commonly implemented using 
digital signal processors. The signal of interest is first processed at the analog front- 
end circuitry, and is then sampled by an analog to digital converter (ADC). 
10 Subsequent processing is implemented on the stored digital samples, which are 
outputted from the ADC. 

The propagation of a signal through a typical wireless channel 
generally results in a number of impairments on the received signal, such as 
multiplicative, dispersive, and additive imperfections. The dispersive imperfections 
15 introduced by multipath propagation may result in inter-symbol interference, which 
can require the use of equalization in the demodulation process. Demodulation refers 
to the extraction of transmitted message data from the samples of the received signal 
passing through the propagation channel. 

Even in the presence of severe impairments, however, algorithms can 
20 be used by the digital signal processor to enable the reliable decoding of the 
transmitted data. 
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One important type of transmission impairment is known as delay 
spread, which is caused by the delay time of a signal being propagated over multiple 
paths (multi-path fading). This delay spread impairment is time-varying in a wireless 
mobile environment and causes interference between differently delayed versions of 
5 the same symbol, which arrive from different paths with varying delay intervals. 

For many practical systems, the decoding to be performed in the case 
of minimal delay spread conditions is relatively simple. On the other hand, the 
decoding process in the presence of larger delay spread conditions requires the use of 
an equalizer scheme. However, the implementation of an equalizer scheme results in 
10 performance variations, and also further complicates the processing required, which 
translates to increases in instruction cycles and power consumption. Therefore, it is 
important to select the most efficient decoding scheme available for the particular 
multi-path channel environment. 

In order to select the most appropriate decoding (demodulation) 
15 scheme, a delay spread estimation method is required. The delay spread estimation 
establishes the essential criterion upon which the selection of a demodulation 
technique is based. 

Known delay spread estimation techniques are generally related to the 
estimation of either the variation of timing and conclusions based on the amount of 
20 this variation, or they involve performing the demodulation under different schemes, 
and concluding that the scheme with the best performance is most appropriate for the 
receive processing function. 

One such scheme is disclosed in US patent 5,400,368. However, the 
disclosed scheme does not take any advantage of the shape and roll-off of the 
25 correlation of the received signal with the synchronization pattern. In addition, the 
disclosed scheme does not take into consideration the effects of different degrees of 
delay spread conditions. 
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Accordingly, it is an object of the present invention to overcome the 
disadvantages of the prior art through an improved method of delay spread estimation, 
which enables the optimal selection of a demodulation technique. 

BRIEF SUMMARY OF THE INVENTION 

In accordance with an illustrative embodiment of the present invention, 
5 a method for estimating the delay spread of a received signal waveform comprises the 
following steps: 

1) receiving a sampled segment of the signal waveform which 
corresponds to a predetermined synchronization sequence; 

2) resampling the received sampled segment at a high sampling rate; 

10 3) correlating the synchronization sequence with the resampled 

segment corresponding to different sampling times; 

4) computing a correlation energy profile from these correlations, and 
determining the peak energy value and its position within the correlation energy 
profile; 

15 5) computing correlation energy values at preselected equal offsets, 

before and after the position of the peak energy value; 

6) determining the maximum of these offset correlation energy values, 
and computing a delay criterion, based on the ratio of this maximum offset correlation 
energy value to the peak energy value; 

20 7) smoothing the data by averaging; 

8) referencing the smoothed delay criterion to a predetermined lookup 

table; and 

9) obtaining the estimated delay spread from the lookup table. 
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In most cases, it is desirable to know if the delay spread falls within a 
particular range. For this purpose, thresholds can be set up in the lookup table to 
identify delay spread ranges. 

An illustrative embodiment of the present invention is more fully 
5 described below in conjunction with the following drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a graph of a typical correlation energy profile. 

Figure 2 is a flow diagram of the inventive delay spread estimation 

method. 

Figure 3 is a flow diagram of a system application of the inventive 

10 method. 

DETAILED DESCRIPTION OF THE INVENTION 

This invention is applicable to digital communication systems that 
employ a synchronization sequence to aid in timing, or initial training. In wireless 
applications, such situations are commonly found in TDMA (Time Division Multiple 
Access) and some CDMA (Code Divisional Multiple Access) systems, since the 
15 timing is critical for the detection of the transmitted data in a mobile environment. 

In one embodiment of the present invention, a synchronization 
sequence has a length Z,, which is denoted by S[i], where i - 0,1,2,...L-1. The receiver 
is designed so that the analog waveform received is sampled at a frequency that 
satisfies or exceeds the Nyquist sampling rate, so as to preclude any aliasing effects. 
20 This implies that the received waveform may be resampled in the digital domain at 
any arbitrary frequency, and with any desired initial sampling phase. 

In one embodiment, resampling of the segment of the received signal 
waveform corresponding to the synchronization sequence is performed at a high 
sampling rate; for example, M times the symbol rate, where Mis chosen in accordance 
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with system implementation requirements. Oversampled received signal samples are 
denoted by r[k], where k = 0,1,2,... The known synchronization pattern is then 
correlated with the resampled versions corresponding to different sampling times, to 
obtain the desired correlations. The computation of a yth correlation result C(j) is 
5 performed in accordance with a relationship written as: 

C(y) = XS*[''M'A* + y] (D 

A peak value of the correlation energy |CG)| 2 is determined, and both a 
location and value of this peak are stored. These parameters are denoted by j pca k and 
|C(j P eak)|\ respectively. In addition, correlation energy values are computed at offsets 
10 of ±A from j pea k, where A is chosen in accordance with a modulation pulse shape, and 
is less than one symbol duration. A representative correlation energy profile is shown 
in Figure 1, with the above described j -values indicated. 

A maximum of the two offset correlation energy values is compared to 
the peak correlation energy value, and the resulting ratio determines the value of the 
1 5 delay criterion D cr ib in accordance with a relationship written as: 

D crit =max(|C(7^ -A)| 2 ,|C( 7 ^ + A)| 2 )/|CO^)| 2 (2) 

Delay criterion value (D c m) can be averaged over multiple received 
synchronization sequences, and increases monotonically with the channel-induced 
delay spread. As such, in one embodiment, this relationship is used in the 
20 determination of delay spread. Since it is usually important to know if the delay 
spread is within a certain range, in some embodiments, suitable thresholds in lookup 
tables are used for decision making. 

A flow diagram of the above described inventive process is shown in 
Figure 2. First, a segment of a signal waveform corresponding to the synchronization 
25 sequence is received 1 . The received segment is resampled 2 at a high sampling rate 
and correlated 3 with a known synchronization pattern. Next, the location and peak 
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value of the correlation energy is found 4 and correlation energy values are obtained 5 
at preselected equal offsets, before and after the peak value. The delay criterion is 
next determined 6, based on a ratio of the maximum offset value and the peak value. 
The delay criterion is averaged 7 over multiple received synchronization sequences 
5 and resultant data is referenced 8 to predetermined thresholds in a lookup table. The 
delay spread estimate is identified 9. 

Computation of the delay criterion D cr ii, as described above, is based 
on values at offsets of ±A. However, in one embodiment, the process is extended to 
encompass a number of intervals that are multiples of A. This extension allows this 
10 embodiment to perform delay spread estimations for a larger range of delay spread 
characteristics. 

The accurate delay spread estimation technique disclosed herein offers 
a specific approach for determining the appropriate method of demodulation. A 
proper selection of demodulation technique becomes critical in wireless systems when 

1 5 optimizing a balance between power consumption and receiver performance. Under 
typical mobile communications applications, it is not uncommon for a mobile 
terminal to experience both high and low delay spread conditions. Therefore, 
satisfying both high and low delay spread conditions, while at the same tome 
efficiently managing power consumption, represents a distinct performance advantage 

20 over the prior art. 

Operation of one embodiment of a system is illustrated by the flow 
diagram of Figure 3. Received sample data is processed 10 as described above to 
determine 20 an accurate delay spread estimate. If a low delay spread is estimated, a 
differential decoding demodulation technique is selected 30. If a high delay spread is 
25 estimated, an MLSE (Maximum Likelihood Sequence Estimation) equalizer 
demodulation technique is selected 40. Finally, an optimally demodulated result is 
output 50 as the decoded result. 



The present invention is applicable to other forms of receiver 
implementation and channel characterization applications, due to its estimation nature 
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and its use of a known characteristic pattern. For example, the present invention can 
be adapted for application to fixed wireless products over certain types of terrain, or to 
any device that performs characterization or testing of a transmission system. 

In short, a method of accurately estimating the delay spread of a multi- 
path propagation channel is disclosed. As a result, an optimal demodulation technique 
is readily determined. Moreover, the disclosed invention can be adapted to a wide 
variety of radio transmission and receiver applications. 

The above described embodiments of the invention are intended to be 
illustrative only. Numerous alternative embodiments may be devised by those skilled 
in the art without departing from the spirit and scope of the following claims. 
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WHAT IS CLAIMED IS: 

1. A method for estimating a delay spread of a received signal 
waveform transmitted through a time-varying multi-path propagation channel, said 
method comprising the steps of: 

sampling a segment of said signal waveform that corresponds to a 
5 predetermined synchronization sequence, 

resampling said sampled received segment at a high sampling rate, 

correlating said synchronization sequence with said resampled segment 
corresponding to different sampling times, 

computing a correlation energy profile from said correlations, 

10 determining the peak energy value and its position within said 

correlation energy profile, 

computing correlation energy values at preselected equal offsets, before 
and after said position of said peak energy value, 

determining a maximum of said offset correlation energy values, 

15 computing a delay criterion, based on the ratio of said maximum offset 

correlation energy value to said peak energy value, and 

referencing said delay criterion to a predetermined lookup table to 
obtain an estimate of the delay spread. 

2. A method in accordance with claim 1 wherein computing a delay 
20 criterion comprises the step of computing a delay criterion averaged over multiple 

received synchronization sequences. 
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3. A method in accordance with claim 1 wherein referencing the delay 
criterion to a predetermined lookup table comprises referencing the delay criterion to a 
predetermined lookup table containing predetermined delay spread thresholds. 

4. A method in accordance with claim 1 wherein computing a 
5 correlation energy profile comprises computing a correlation energy profile written as: 

C(y) = 2S*[/MiM + 7] 

where C(j) is a yth correlation result; 

S[i] is a synchronization sequence, where /=0,1,2,...L-1; ' 

L is a length of the synchronization sequence; 

10 r[k] is a sequence of oversampled received signal samples, where 

A=0,l,2,...; and 

a resampling rate of the received signals is M times a symbol rate of the 
received signal. 

5. A method in accordance with claim 1 wherein computing a delay 
1 5 criterion comprises computing a delay criterion written as: 

D cnl =max(|C( 7 ^ -A)| 2 ,|CO^ + A)| 2 )/|C(;^)| 2 
where D C rit is the delay criterion; 

Cijpeak) is a peak value of correlation energy |C(/)| 2 , where C{j) is a yth 
correlation result; 

20 jpeak is a location of the peak value of correlation energy; 

and A is an offset less than one symbol duration from selected in 
accordance with a modulation pulse shape. 
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6. A method in accordance with claim 1 wherein computing a delay 
criterion comprises the step of computing a delay criterion in accordance with 
multiples of said offsets. 

7. A method in accordance with claim 1 further comprising the step of 
5 selecting a demodulation technique in accordance with the estimated delay spread. 

8. A method in accordance with claim 7 wherein selecting a 
demodulation technique in accordance with the estimated delay spread comprises the 
step of selecting a selecting a differential decoding demodulation technique when 
estimated delay spread is below a predetermined threshold. 

10 9. A method in accordance with claim 7 wherein selecting a 

demodulation technique in accordance with the estimated delay spread comprises the 
step of selecting an MLSE (Maximum Likelihood Sequence Estimation) equalizer 
demodulation technique when the estimated delay spread is above a predetermined 
threshold. 

15 10. A communication receiver device having delay spread estimation 

capability, comprising: 

a receiver circuit for receiving a sampled segment of a signal waveform 
which corresponds to a predetermined synchronization sequence, 

a resampling circuit for digitally sampling said received sampled 
20 segment at a high sampling rate, 

a processing circuit for correlating said synchronization sequence with 
said resampled segment corresponding to different sampling times, 

said processing circuit for computing a correlation energy profile from 
said correlations, 

25 said processing circuit for determining the peak energy value and its 

position within said correlation energy profile, 
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said processing circuit for computing correlation energy values at 
preselected equal offsets, before and after said position of said peak energy value, 

said processing circuit for determining a maximum of said offset 
correlation energy values, 

5 said processing circuit for computing a delay criterion, based on the 

ratio of said maximum offset correlation energy value to said peak energy value, and 

a comparison circuit for determining said delay spread estimation by 
referencing said delay criterion to a lookup table. 

11. The communication device of claim 10 wherein said delay 
10 criterion is averaged over multiple received synchronization sequences. 

12. The communication device of claim 10 wherein said computation 
of said correlation energy profile in step d) is derived from an equation written as: 

where C(j) is a yth correlation result; 
1 5 S[i] is a synchronization sequence, where /=0,1 ,2,...Z,-1 ; 

I is a length of the synchronization sequence; 

r[k] is a sequence of oversampled received signal samples, where 

A=0,l,2,...;and 

a resampling rate of the received signals is M times a symbol rate of the 
20 received signal. 

13. The communication device of claim 10 wherein said computation 
of said delay criterion in step h) is derived from an equation written as: 
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D cru = max(\C(j peak -A)|\|C(y^ +A)| 2 )/|C(y^)| 2 
where £> C rit is the delay criterion; 

CQ'peak) is a peak value of correlation energy |C(/)| 2 , where C(j) is a jth 
correlation result; 

jpeak is a location of the peak value of correlation energy; 

and A is an offset less than one symbol duration from jpeak selected in 
accordance with a modulation pulse shape. 

14. The communication device of claim 10 wherein said computation 
of said delay criterion is based on multiples of said offsets. 
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